Dissipative Macroscopic Quantum Tunneling in Type-I Superconductors 
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We study macroscopic quantum tunneling of interfaces separating normal and superconducting 
regions in type-I superconductors. Mathematical model is developed, that describes dissipative 
quantum escape of a two-dimensional manifold from a planar potential well. It corresponds to, 
e.g., a current-driven quantum depinning of the interface from a grain boundary or from artificially 
manufactured pinning layer. Effective action is derived and instantons of the equations of motion 
are investigated. Crossover between thermal activation and quantum tunneling is studied and the 
crossover temperature is computed. Our results, together with recent observation of non-thermal 
low-temperature magnetic relaxation in lead, suggest possibility of a controlled measurement of 
quantum depinning of the interface in a type-I superconductor. 

PACS numbers: 74.25.Ha, 74.50.+r, 75.45.+j, 03.75.Lm 



I. INTRODUCTION 

Macroscopic quantum tunneling refers to the situa- 
tion when an object consisting of many degrees of free- 
dom, coupled to a dissipative environment, escapes from 
a metastable well via underbarrier quantum tunneling 1 . 
In condensed matter this phenomenon was first observed 
through measurements of tunneling of the macroscopic 
magnetic flux created by a superconducting current in a 
circuit interrupted by a Josephson junction 2 . Another 
example is tunneling of magnetization in solids 3 . In 
cases of the magnetic flux or the magnetic moment of a 
nanoparticle, the tunneling object is described by one or 
two macroscopic coordinates that depend on time, like 
in a problem of a tunneling particle in quantum me- 
chanics. The environment enters the problem through 
interaction of these macroscopic coordinates with micro- 
scopic excitations of the medium. Equally interesting, 
but significantly more involved, is the problem of tunnel- 
ing of a macroscopic field between two distinct configu- 
rations. Examples are tunneling of vortex lines in type- 
II superconductors^^ and tunneling of domain walls in 
magnets^—. The essential difference between the last two 
examples is that tunneling of vortex lines is determined 
by their predominantly dissipative dynamics^— , while 
tunneling of the spin-field is affected by dissipation to 
a much lesser degree. Theory that describes quantum 
tunneling of extended condensed-matter objects involves 
space-time instantons that are similar to the instantons 
studied in relativistic field models. Examples that are 
available for experimental studies are limited. Conse- 
quently, any new example of tunneling of an extended 
object must be of significant interest. 

Recent measurements of low-temperature magnetic 
relaxation of leadi 5 . have elucidated the possibility of 
macroscopic quantum tunneling in type-I superconduc- 
tors. Such superconductors (with lead being a prototyp- 
ical system), unlike type-II superconductors, do not de- 



velop vortex lines when placed in the magnetic field. In- 
stead, they exhibit intermediate state in which the sam- 
ple splits into normal and superconducting regions sepa- 
rated by planar interfaces of positive energy^—. Equi- 
librium states and dynamics of interfaces have been well 
studied by now-i^r— . In all these studies the interface was 
treated as a classical object. Recently, however, it was 
noticecU 5 - that slow temporal evolution of magnetization 
in a superconducting Pb sample was independent of tem- 
perature below a few kelvin. This observation pointed 
towards possibility of quantum tunneling of interfaces in 
the potential landscape determined by pinning. In gen- 
eral the pinning potential would be due to random dis- 
tribution of pinning centers or due to properties of the 
sample surface. In a polycrystalline sample it may also 
be due to extended pinning of interfaces by grain bound- 



Modern atomic deposition techniques permit prepara- 
tion of a pinning layer with controlled properties. This 
inspired us to study a well defined problem in which 
the interface separating normal and superconducting re- 
gions is pinned by a planar defect. The corresponding 
pinning barrier can be controlled by a superconduct- 
ing current that exerts a force on the interface. At low 
temperature the depinning of the interface would occur 
through quantum nucleation of a critical bump shown 
in Fig. [TJ Somewhat similar problems in 1+1 dimen- 
sions have been studied for a flux line pinned by the in- 
terlayer atomic potential in a layered superconductor— 
and for a flux line pinned by a columnar defect 25 . How- 
ever, the two-dimensional nature of the interface, as com- 
pared to a one-dimensional flux line, makes the inter- 
face problem more challenging. Note that tunneling of 
two-dimensional objects has been studied theoretically in 
application to non-thermal dynamics of planar domain 
walls 7 and quantum nucleation of magnetic bubbles 2 ^. 
These studies employed non-dissipative dynamics of the 
magnetization field because corrections coming from dis- 
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of the interface associated with its total area, the sec- 
ond term is due to the space-dependent potential energy, 
V [x, y, Z(x, y)] , of the interface inside the imperfect crys- 
tal, and the third term is due to dissipation 1 . Same as for 
the flux lines, we neglect the inertial mass of the interface. 
Its dynamics in a type-I superconductor is dominated by 
friction. 

We consider pinning of the interface by a planar defect 
located in the XY plane and choose the pinning potential 
in the form 



Figure 1: Interface between normal and superconducting re- 
gions in a type-I superconductor, pinned by a planar defect 
in the XY plane. Transport current parallel to the interface 
controls the energy barrier. Depinning of the interface occurs 
through quantum nucleation of a critical bump described by 
the instanton of the equations of motion in 2+1 dimensions. 



sipation are not dominant for spin systems. On the 
contrary, the Euclidean dynamics of the interface in a 
type-I superconductor is entirely dissipative, described 
by integro-differential equations in 2+1 dimensions. As 
far as we know this problem has not been studied before. 

The article is structured as follows. Theoretical model 
is formulated in Sec. [TXJ Properties of the pinning poten- 
tial and the effective action in the vicinity of the critical 
depinning current are analyzed in Sec. Mil Instantons 
of the dissipative model in 2+1 dimensions are investi- 
gated in Sec. IIVI Crossover from quantum tunneling to 
thermal activation is studied in Sec. [V] Sec. I VII contains 
estimates of the effect and final conclusions. 



II. THE MODEL 

We describe the interface by a smooth function Z(x, y), 
see Fig. [TJ Dimensionless Euclidean effective action as- 
sociated with the interface is 
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pal dxdy (--r-^ 
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where 2a is roughly the width of the well that traps the 
interface and p < 1 is a dimensionless constant describing 
the strength of the pinning. The interface separates the 
normal state at Z < from a superconducting state at 
Z > 0. Superconducting current parallel to the planar 
defect (and to the interface pinned by the defect) exerts a 
Lorentz force on the interface similar to the force acting 
on a vortex line in a type-II superconductor. We shall as- 
sume that the magnetic field is applied in the y direction 
and that the transport current of density j flows in the 
x direction. The driving force experienced by the dxdy 
element of the interface in the z direction is given by 



d 2 F z 
dxdy 



= - c J dzjB(z) 



(4) 



Here B(z) — B c exp(—z/5) is the magnetic field inside 
the interface with S = y/^X£. Integration then gives 
d 2 F z j (dxdy) — B c Sj/c. The corresponding contribu- 
tion to the potential can be obtained by writing F z as 
-VzVr,, yielding 



d 2 V L (Z) B c 8 



dxdy 

The total potential, V(Z) 



jZ. 



V P (Z) + V L {Z) is 



(5) 



Seff = j. f dr / dxdy [l + (VZ) 2 ]" 



dr / dxdyV [x,y,Z(x,y,T)} 



(1) 



i dr / dr' / dxdy 



[Z(x,y,r) - Z(x,y,T')Y 



r /\2 



where r is the imaginary time, a is the surface energy 
density of the interface and 77 is a drag coefficient, given 
respectively b y 15 ' 27 



V(Z) = pa / dx dy 



-J*+t-t) < 6) 



where we have introduced dimensionless Z = Z/a and 
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with k = Al/£. Note that for a type-I superconductor 
k < 1/V2. 



V = 



2p n c 2 



(2) 



with B c being the thermodynamic critical field, £ be- 
ing the superconducting coherence length, Al being the 
London length, and p n being the normal state resistiv- 
ity. The first term in Eq. {1} is due to the elastic energy 



III. EFFECTIVE ACTION IN THE VICINITY 
OF THE CRITICAL CURRENT 

Measurable quantum depinning of the interface can oc- 
cur only when the transport current is close to the crit- 
ical current, j c , that destroys the energy barrier. It is, 
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therefore, makes sense to study the problem at j — s- j c 
Maxima and minima of the function 



f(lZ) 



-]Z 
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(8) 



that enters Eq. (|6j) are given by the roots of the equa- 
tion Z 3 - Z + j = 0. At j 2 < 4/27 it has three real 
roots corresponding to one minimum and two maxima of 
the potential on two sides of the pinning layer, whereas 
at j > 4/27 there is one real root corresponing to the 
maximum of /. Consequently, the barrier disappears at 
j 2 = 4/27, providing the value of the critical current 



Jc 



3V3' 



3c 



2pcB c 



97r\/6/ia 



(9) 



At j = j c the minimum and the maximum of the po- 
tential combine into the inflection point Z = Z c given by 
the set of equations 



= -Zi 







-3Zi 



Zc - jc 

f 1 



(10) 



that correspond to zero first and second derivatives of /. 
The value of Z c deduced from these equations is l/\/3- 
It is convenient to introduce small parameter 



e = 1 - j/jc 
so that j = j c (l — e) and 

- - , ^ 2 
j = jc(l ~ e) 



3V3 



(1-e) 



(11) 



(12) 



Let Zq(J) be the minimum of / (see Fig. [5]) satisfying 

Z 3 - Z Q +] c (l-e) = 0. (13) 

Consider Z' = Z — Zq. It is easy to find that the form of 
the potential in the vicinity of Zq is 



f = .f[Z Q)] + )-{!- 3Z 2 )Z 



2* 



(14) 



At small e one has Zq — > Z c = l/\/3, so that 1 — 3Z'q in 
front of Z' 2 in Eq. (fT4"| is small. The first term in Eq. 
(TH)) can be omitted as unessential shift of energy, while 
the last term proportional to Z' 4 can be neglected due 
to its smallness compared to other Z'-dependent terms. 
Consequently, one obtains the "effective potential" 

feffti, Z) = ~(1 - 3Z 2 )Z 2 - Z Z 3 . (15) 

We need to know the dependence of Zq on e. Writing 
Z (e) = Z c [l-/3(e)], with the help of Eq. CE3]), we obtain 
/3(e) = y/2e/3 to the lowest order on e. Then 1 - 3Z$ « 
2y/2e/3 and 



feff(e,Z) 



Z*_ 

V3 



(16) 




Figure 2: Effective potential 

The height of the effective potential is "y^ 6 "^ 2 anc l * ne 

width is \/2e, see Fig. [2] 

As follows from the equations of motion, smallness of 
e results in [VZ\ ~ pe <C 1. This allows one to replace 

[l + (VZ) 2 ] 5 in Eq. (J) with 1 + \{V Z) 2 . Introducing 
dimensionless variables 

(2 v JitB 2 c \_ f _ ( / _ xi/a(a. >y ) 



.r 



V 3V3 m 2 

v = V(x,y,Z)/crp, u 



t, (x 1 ,x 2 ) = (V 2e / 3 p) 



Z/a - Z c (l - V2e/3) 



we obtain 



S, 



eff 



ye ya 
3\/6iTp fi 



dxo / dx± dx 2 



(17) 



1 { , [u(x > xi,x 2 ) - u(x' ,xi,x 2 )] 

O.Xrt 



R 



(x - x' ) 2 



(18) 



where V = (8i,8 2 ). 



IV. INSTANTONS OF THE DISSIPATIVE 2+1 
MODEL 

Quantum depinning of the interface is given by the in- 
stanton solution of the Euler-Lagrange equations of mo- 
tion of the 2+1 field theory described by Eq. (fTgj) : 



y — 



SC 



5 (8u/8x^) 



8u 



(19) 



This gives 



V^u-2u+u -2 / dx' 



, u(x ,xi,x 2 ) - u(x' a ,xi,x 2 ) 



(xq - x' ) 2 



= 
(20) 
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with the boundary conditions 

u(— Q/2,xi,£ 2 ) = u(Q/2,xi,x 2 ) V(xi,x 2 ) € R 2 

max u(xo, x\, x 2 ) = u(0, x%, x%) V(ii,i2) eR 2 

xoe[— f2/2,Q/2] 

(21) 

that must be periodic on imaginary time with the period 
h/(ksT). The corresponding period on xq is 



The output is the pair (A„, u n (uj)). Finally, we apply 
a rescaling of u n by a factor (27r) 3 / 2 A„ to obtain the 
instanton solution. This procedure leads to 



b eff = —+0 



(2 P yTeZB 2 c \ ft 
V 3^3 W 2 J k B T ' 



(22) 



This equation cannot be solved analytically, so we must 
proceed by means of numerical methods. 



A. Zero temperature 



We apply the Fourier transform 
1 



u(uj) 



(2tt) 3 / 2 
to equation (|2T)|) and get 

(2tt)- 3 / 2 



u(x)e 



(23) 



u(uj) 



2 + 27t|cj | + u>{ + uj 2 7r3 



d oj'u{w - i3')u(C3') 
(24) 

which is still an integral equation for u(u). The effective 
action JT8]) in terms of u(qj) becomes 



eff M 



3^/6np ft 



R 3 



d 3 uj u(uj)u(—uj) x 



+ wf) + 1 + tt\uj q \ 



3(2tt) 3 / 2 



d 3 w & a uj u(u)u(uj )u(-uj - a; ) 



.(25) 



We use the algorithm that is a field-theory extension of 
the algorithm introduced in Refs. |28||29| for the problem 
of dissipative quantum tunneling of a particle. It consists 
of the following steps: 

1. Start with an initial aproximation uq((2). Define 
the operator 



Sy/Gnp ft 



(27) 



with numerical value of the integral Iq — 531 ± 19. This 
somewhat surprisingly large value of the integral has been 
confirmed by our use of different computational grids. 



B. Non-zero temperature 

At T 7^ the period of the instanton solution is finite, 
given by Eq. (|2"2"|) . We look for a solution of the type 



u(x 



,X!,X 2 ) = ^ ' 



>u n (xi,x 2 ) 



(28) 



nez 



with LUg.n — 2nn/Q. Introducing into (|20p the above 
functional dependence and applying a 2D Fourier trans- 
form we obtain 



u„(uj) 



1 



2 + 27t|cj 0i „| + uj 2 
r~y^ / d 2 u/w„_ p (u; - ti')u p {w') ) , 

i7r r^L J-r? I 



(29) 



which is the integral equation for u n with uj = (uji,uj 2 ). 
In terms of {u n (u;)} the effective action becomes 



(30) 



V" / d 2 uju n (uj)u- n (-Lj) I — + 1 + 7r|cj ,«| 
— ^ / d 2 wd 2 w'ii„(w)ii ra (w')ii-„- m (-w - w') 



()7T 



ti.m^Z 



The numerical algorithm is analogous to the one used in 
the T = case. It leads to 



5', 



0:Rx£ 2 (R 3 

(A,u(w)) i— > 



£ 2 (R 3 ) 
A 



(26) 



3\/67rp ft 



J(T) 



(31) 



3 / ^ The value of the integral depends on the value of T in 

d uj u(ui — uj )u(uj ) corn p ar i son w ith the temperature, T c , of the crossover 
from quantum tunneling to thermal activation (see be- 
low). At T <C T c the numerical value of I(T) is very 
close to Iq, while at T > T c we recover the Boltzmann 
exponent, S e ff = Vo/QcbT), with Vb being the energy 
barrier for depinning. Computation of I(T) in the inter- 
mediate temperature range requires very large computer 
time and will be reported elsewhere. Nevertheless, as 
5. Repeat steps (2) — (4) until the successive difference we shall see below, the crossover temperature T c can be 



2 + 2it\uj \ +uj\+ujI J r3 

2. Let U\{uj) = 0(Xq,Uo(uj)) for an initial A € R 

3. Calculate Ai = A /£ 2 with £ 

4. Find u 2 (uj) = 0(Xi,ui(uj)). 



ui(uj— 0) 

«o(<3=0) ' 



satisfies a preset convergence criterion. 



computed exactly. 
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CROSSOVER TEMPERATURE 



VI. DISCUSSION 



The crossover temperature can be computed by means 
of theory of phase transitions^. Above T c , the so- 
lution minimizing the instantion action is a function 
u(xq, xi, X2) = Uo(xx,X2) that does not depend on ad- 
just below T c , the instanton solution can be split into 
the sum of uq and a term that depends xq, 



u(x ,xi,x 2 ) = u (xi,x 2 ) + ui(xi,X2)cos(27r/f2x ) 



The instanton action is proportional to 

/ dxi dx2&(xi, X2',u,Vu) . 

Js. 2 



(32) 



(33) 



where <&(xi , a^i u, Vu) is the spatial action density. Using 
the expansion of u introduced in the previous section, we 
obtain 



n 



n 



(Vu ) 2 + v(u ) 



-(v Ul y + Aui + 0(4) 

with v(it) = « 2 — u 3 /3 and 

A=^"(U )+7T 2 



(34) 



(35) 



If A > 0, the only (u±, Vui) minimizing $ is iti = 0, 
so we define the crossover temperature by the equation 



min A = min — -v"\un(xi . X2)] + 7r 2 
sm 2 xgr 2 4 



0. 



(36) 



Notice that this minimum corresponds to the minimum 
of v"[uo(xi, X2))]. The equation of motion for uq is 



V 2 u - 2u + Gq = . 



(37) 



Solution corresponding to the minimum is spherically 
symmetric, 



u = u \r 



(38) 



satisfying boundary conditions: uq — > at r — > 00 and 
Mo(0) = 3, which is the width of the potential. Conse- 
quently, 



min v"[u (xi, x 2 )1 = min v"(u ) 
fei 2 a e[o,3] 

= min 2(1 - U ) = -4 . 

« e[o,3] 



(39) 



Then, according to equations (|55|) and (pJB")) , the crossover 
temperature is determined by the equation f2(T c ) = 7r 2 , 



which gives 



3VSn 2 k B ria 2 3tt 2 V3k k B a? 



We are now in a position to discuss feasibility of the 
proposed experiment on quantum dcpinning of the inter- 
face from a planar defect in a type-I superconductor. Two 
conditions must be satisfied. Firstly the dimensionless ef- 
fective action of Eq. (f27|) . which is the WKB exponent of 
the tunneling rate, should not exceed 30 — 40 in order for 
the tunneling to occur on a reasonable time scale. Sec- 
ondly, the crossover temperature determined by Eq. (|40[) 
better be not much less than one kelvin. For a known 
superconductor, the two equations contain three param- 
eters: The parameter p < 1 describing the strength of 
pinning, the parameter a describing the width of the pin- 
ning layer, and the parameter e that controls how close 
the transport current should be to the depinning cur- 
rent. We, therefore, have to investigate how practical is 
the range of values of these parameters that can provide 
conditions S e ff ~ 30 and T c ~ IK. 

Let us choose lead as an example. The values of Al 
and £ in lead are 37 nm and 38 nm, respectively, giving 
« = = 0.45. The critical field is B c w 800 G. The 

elastic energy of the interface is a « 0.4 erg/cm 2 . The 
normal state resistivity in the kelvin range is 5 x 10~ u 
f2-m = 5.6 x 10 _21 s., while the drag coefficient is 77 w 
0.35 erg-s/cm 4 . Then equations (|27|) and (|40[) with con- 
ditions S e ff ~ 30 and T c ~ IK give a/p 1 ^ 3 ~ 3.7 nm 
and y/ea ~ 0.25 nm. If the pinning layer is incompatible 
with superconductivity, then at 2a < £ one should expect 
p ~ 2a/£, giving a ~ 1.65 nm and e ~ 0.02. This means 
that observation of quantum escape of the interface from 
a pinning layer of thickness 2a ~ 3.3 nm in a supercon- 
ducting Pb sample at T ~ IK would require control of 
the transport current within two percent of the critical 
depinning current. All the above parameters are within 
experimental reach. 

In conclusion, we have studied quantum escape from a 
planar pinning defect of the interface separating super- 
conducting and normal regions in a type-I superconduc- 
tor. This can correspond to either quantum depinning of 
the interface from a grain boundary or quantum depin- 
ning from an artificially prepared layer. The computed 
tunneling rate, the required temperature and other pa- 
rameters all fall within realistic experimental range. We 
encourage such experiment as it would present a rare op- 
portunity to study, in a controllable manner, dissipative 
quantum tunneling of an extended object. 



VII. ACKNOWELEDGEMENTS 

The work at the University of Barcelona has been 
supported by the Spanish Government project No. 
MAT2008-04535 and by Catalan ICREA Academia. R.Z. 
acknowledges financial support from the Ministerio de 
Ciencia e Innovation de Espana. The work of E.M.C. 
at Lehman College has been supported by the U.S. 
Department of Energy through grant No. DE-FG02- 



6 



93ER45487. 



1 A. O. Caldeira and A. J. Leggett, Phys. Rev. Lett. 46, 211 
(1981); Ann. Phys. (N.Y.) 149, 374 (1983). 

2 J. Clarke, A. N. Cleland, M. H. Devoret, D. Esteve, and J. 
M. Martinis, Science 239, 992 (1988). 

3 E. M. Chudnovsky and J. Tejada, Macroscopic Quantum 
Tunneling of the Magnetic Moment (Cambridge University 
Press, Cambridge, England, 1998). 

4 G. Blatter, M. V. Feigel'man, V. B. Geshkenbein, A. I. 
Larkin, and V. M. Vinokur, Rev. Mod. Phys. 66, 1125 
(1994). 

5 Y. Yeshurun, A. P. Malozemoff, A. Shaulov, Rev. Mod. 
Phys. 68, 911 (1996). 

6 T. Nattermann and S. Scheidl, Adv. Phys. 49, 607 (2000). 

7 E. M. Chudnovsky, O. Iglesias, and P. C. E. Stamp, Phys. 
Rev. B 46, 5392 (1992). 

8 K. Hong and N. Giordano, J. Phys.: Cond. Matter 8, L301 
(1996). 

9 J. Brooke, T. F. Rosenbaum, and G. Aeppli, Nature 413, 
610 (2001). 

10 G. Blatter, V. B. Geshkenbein, and V. M. Vinokur, Phys. 
Rev. Lett. 66, 3297 (1991). 

11 B. I. Ivlev, Yu. M. Ovchinnikov, R. S. Thompson, Phys. 
Rev. B 44, 7023 (1991). 

12 J. Tejada, E. M. Chudnovsky, and A. Garcia, Phys. Rev. 
B 47, 11552 (1993). 

13 P. Ao and D. J. Thouless, Phys. Rev. Lett. 72, 132 (1994). 

14 M. J. Stephen, Phys. Rev. Lett. 72, 1534 (1994). 

15 E. M. Chudnovsky, S. Velez, A. Garci'a-Santiago, J. M. 
Hernandez, and J. Tejada, Phys. Rev. B 83, 064507 (2011). 

16 L. D. Landau, Phys. Z. Sowjetunion 11, 129 (1937). 



17 Y. V. Sharvin, Zh. Eksp. Teor. Fiz. 33, 1341 (1957). 

18 R. P. Huebener, Magnetic Flux Structures of Superconduc- 
tors ( Springer- Verlag, New York, 1990). 

19 A. V. Kuznetsov, D. V. Eremenko, and V. N. Trofimov, 
Phys. Rev. B 57, 5412 (1998). 

20 A. Cebers, C. Gourdon, V. Jeudy, and T. Okada, Phys. 
Rev. B 72, 014513 (2005). 

21 M. Menghini and R. J. Wijngaarden, Phys. Rev. B 72, 
172503 (2005). 

22 R. Prozorov, Phys. Rev. Lett. 98, 257001 (2007). 

23 R. Prozorov, A. F. Fidler, J. R. Hoberg, and P. C. Canfield, 
Nature Physics 4, 327 (2008). 

24 S. Velez, A. Garcfa-Santiago, J. M. Hernandez, and J. Te- 
jada, Phys. Rev. B 80, 144502 (2009). 

25 E. M. Chudnovsky, A. Ferrera, and A. Vilenkin, Phys. Rev. 
B 51, 1181 (1995). 

26 E. M. Chudnovsky and L. Gunther, Phys. Rev. B 37, 9455 
(1988); A. Ferrera and E. M. Chudnovsky, Phys. Rev. B 
53, 354 (1996). 

27 E. M. Lilshitz and L. P. Pitaevskii, Statistical Physics, Part 
2 (Oxford: Pergamon, 1980). 

28 L. Chang and S. Chakravarty, Phys. Rev. B 29, 130 (1984). 

29 D. Waxman and A. J. Leggett, Phys. Rev. B 32, 4450 
(1985). 

30 I. Affleck, Phys. Rev. Lett. 46, 388 (1981); A. I. Larkin and 
Yu. N. Ovchinnikov, Pis'ma Zh. Eksp. Teor. Fiz. 37, 322 
(1983) [JETP 37, 382 (1983)]; E. M. Chudnovsky, Phys. 
Rev. A 46, 8011 (1992). 



